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Abstract 

We calculate at two-loop order the complex-valued scattering amplitude related to 
the twice-iterated scalar-isovector boson-exchange between nucleons. In comparison to 
the once-iterated boson-exchange amplitude it shows less dependence on the scattering 
angle. We calculate also the iteration of the (static) irreducible one-loop potential with 
the one-boson exchange and find similar features. Together with the irreducible three- 
boson exchange potentials and the two-boson exchange potentials with vertex corrections, 
which are also evaluated analytically, our results comprise all nonrelativistic contributions 
from scalar-isovector boson-exchange at one- and two-loop order. The applied methods 
can be straightforwardly adopted to the pseudoscalar pion with its spin- and momentum- 
dependent couplings to the nucleon. 
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Over the last years, effective field theory (in particular chiral perturbation theory) has been 
successfully applied to the two- nucleon system at low and intermediate energies [1, 2, 3]. The 
strategy is to construct the long- and medium-range NN-potential systematically from one- 
, two- and three-pion exchanges and to represent the short- distance dynamics by adjustable 
contact interactions. In order to obey unitarity of the S-matrix the NN-potential must then 
be iterated to all orders, e.g., by solving the Lippmann-Schwinger integral equation. In prac- 
tice this unitarization procedure introduces additional cut-offs or off-shell form factors in order 
to eliminate high-momentum components from the chiral NN-potential. However, there exist 
also purely algebraic methods for unitarization, such as the inverse amplitude expansion or the 
N/D-method, which do not require such additional regularizations. In the context of coupled 
channels dynamics the so-called chiral unitary approach has been widely used for meson-meson 
and meson-baryon scattering [4]. In a recent work [5], Oiler has adopted the method to elastic 
nucleon-nucleon scattering, considering so far only contact interactions. Clearly, a more realistic 
treatment must include the one-pion exchange as well as pion-loop contributions. If one goes 
now in this novel unitarization scheme for NN-scattering to sufficiently high orders in the small 
momentum expansion one will encounter the (multiply) iterated pion-exchanges. At present, 
analytical expressions are only known for the once-iterated pion-exchange scattering amplitude 
(see sect. 4.3 in ref. [6]), but not for the higher iterations or the iterations involving the pion-loop 
potential. 

The purpose of the present paper is to perform such calculations for the technically simpler 
case of a scalar-isovector boson (treating the nucleons in the nonrelativistic approximation). The 
generalization to the pseudoscalar pion with its spin- and momentum-dependent couplings to the 
nucleon is in principle straightforward and will be presented elsewhere [7] together with results 
for NN-phase shifts and mixing angles. For the following we assume that the scalar-isovector 
boson of mass m is light (m <C M, where M denotes the nucleon mass) and that it couples weakly 
to the nucleon (with a coupling constant g ~ 27rm/M). Let us start with recalling the once- 
iterated meson-exchange amplitude. The corresponding planar one-loop diagram for the elastic 
scattering process in the center-of-mass frame is shown in Fig. 1. The pertinent complex- valued 



Figure 1: Diagrams of once-iterated (left) and twice- iterated (right) boson-exchange 



one-loop integral reads: 
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and it can be solved in terms of inverse trigonometric and logarithmic functions 
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Here, p =\p\ = \p'\ is the center-of-mass momentum and q = \p' — p \ denotes the momentum 
transfer between the two nucleons. ri,2 are the usual isospin operators. For the value at threshold 
p = = one easily deduces from Eq.(2): G{0, 0) = l/(2m^). 

In Fig. 2 we show the real and imaginary part of the (dimensionless) once-iterated scalar- 
boson exchange amplitude m^G{p, q) as a function of the center-of-mass momentum p in the 
region < p < 2m. The three full and three dashed lines correspond to the cases of forward 
scattering [z = 1), perpendicular scattering {z = 0), and backward scattering {z = —1) with 



the momentum transfer given hy q — p^2{l — z). All other possible kinematical configurations 
lie of course in between. One observes from Fig. 2 a rapid decrease of the real part, whereas the 
imaginary part develops a broad maximum around p ~ 0.5m. 

Next, we come to the twice-iterated meson-exchange represented by the right diagram in 
Fig. 1. The corresponding two-loop integral: 
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defines the complex- valued scattering amplitude H{p,q). In order to evaluate it, we make use 
of (perturbative) unitarity in the form of the Cutkosky cutting rule. For the problem at hand 
it states that the imaginary part of the twice-iterated boson-exchange amphtude is equal to 
the two-body phase space integral (i.e. a solid angle integral) over the one-boson exchange 
amplitude times twice the real part of the once-iterated boson-exchange amplitude. Since the 




Figure 2: Real and imaginary part of the once-iterated scalar-boson exchange NN-scattering 
amplitude G{p,q) versus the center-of-mass momentum p. The cases z = 1,0,-1 correspond 
to forward, perpendicular, and backward scattering with momentum transfer q = p^j2{l — z), 
respectively. 



latter is known in closed form (see Eq.(2)) we get after some transformations the following 
single-integral representation for the imaginary part: 
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The associated real part is obtained from an unsubtracted dispersion relation of the form: 

Re H{p, q)^-r Im H{p', q) . 
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It is astonishing that the value at threshold p = g = is still calculable analytically, with the 
result: 

m^H{Q, 0) = In ^ . (6) 



In case of the two-loop representation Eq.(3) we can show some intermediate steps: 
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Figure 3: Real and imaginary part of the twice- iterated scalar-boson exchange NN-scattering 
amplitude H{p, q). For further notations, see Fig. 2. 
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Figure 4: Crossed two-boson exchange diagram, and diagrams related (partly) to the iteration 
of the irreducible one- loop potential V2{q) with one-boson exchange. 



which lead to this results. For the numerical evaluation of the principal value integral in Eq.(5) 
it is advantageous to convert it into a sum of nonsingular integrals by the identity: 
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where s = p^. In Fig. 3 we show the real and imaginary part of the (dimensionless) twice- 
iterated scalar-boson exchange amplitude m'^H{p, q) versus the center-of-mass momentum p in 
the interval < p < 2m. In comparison to the once-iterated boson-exchange amplitude G{p, q) 
one observes less dependence on the scattering angle, a feature which is exhibited by the weaker 
splitting of the curves {z = 1) and {z = —1) corresponding to forward and backward scattering. 
Another remarkable property is that the real part Re H(j), q) crosses zero at p ~ 0.58 m and from 
there on it continues with small negative values which asymptotically tend to zero. At two-loop 
order there are also the non-planar ladder diagrams shown in Fig. 4. These correspond (partly) 
to the iteration of the irreducible two-boson exchange (i.e. the one-loop potential) with the 
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one-boson exchange. The complete one-loop potential arises from the (left) crossed two-boson 
exchange diagram in Fig. 4 and the irreducible part of the (left) box diagram in Fig. 1. Both 
pieces together lead to the following result for the (static) one-loop potential^ in momentum 
space: 

V^W = — n •r2^==^ln . 9 

Our sign convention is chosen here such that the (static) one-boson exchange reads: Vi{q) — 
9^ 7^1 ■ '^2 {iTi'^ + Q'^)~^- It is interesting to note here that modulo their isospin factors 3 ± 2ri ■ T2 
the crossed box diagram and the irreducible part of the planar box diagram are equal but with 
opposite sign in the limit M — > 00 (for further details on that fact, see sect. 4.2 in ref.[6]). The 
coordinate space potential corresponding to Eq. (9) can be expressed through a modified Bessel 
function: V2{r) — —g^Ti ■ T2 -f^o(2mr)/27rV. 

Now we are in the position to calculate also the iteration of the (static) one-loop two-boson 
exchange potential with the one-boson exchange: 

n (P, ^) = - 2n • T2) J(p, q) . (10) 

Unitarity determines again the imaginary part of the complex-valued amplitude J(p, q) as: 
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and the real part Re J(p, q) follows from an unsubtractcd dispersion relation analogous to Eq.(5). 
The value at threshold p = g = is again calculable analytically, with the result: 

m3 J(0, 0) = r = vrf 1 - ^ V (12) 

where we have made use of the spectral function representation [6] of the (static) one-loop 
potential V2{q) written in Eq.(9). In Fig. 5 we show the real and imaginary part of the (dimen- 
sionless) iterated one-loop amplitude m^J{p, q) versus the ccntcr-of-mass momentum p in the 
region < p < 2m. The real part displays only a minor dependence on the scattering angle, 
which is more pronounced for the imaginary part. Another noticeable feature is that the real 
part stays positive throughout and that the decrease of the amplitudes with increasing momen- 
tum p is slower in comparison to G{p, q) and H{p, q) shown Figs. 2,3. The latter property has 
to do with the presence of the larger mass scale 2m in the (static) one- loop potential V2(g). At 
two-loop order there are in addition the three crossed thrcc-boson exchange diagrams shown in 
Fig. 6. These build up together with the irreducible parts of the two-loop ladder diagrams in 
Figs. 1,4 the (irreducible) three-boson exchange potential. The separation of a two- loop ladder 
diagram into twice-iterated (~ M^), once-iterated (~ M^), and nonrelativistic irreducible com- 
ponents (~ Al^) is defined by the expansion of the corresponding Feynman amplitude in powers 
of the reciprocal nuclcon mass 1/Af. For that one integrates first the product of light boson 
propagators (/| — rn? + iO+)~^ and heavy nucleon propagators 2M{p'^ — + iO^)^^ over the 
loop-energies Z°2 via residue calculus and then expands in powers of 1/M. This way we find 



^In this work, we do not consider relativistic 1/M-corrections to the irreducible one-, two- and three-boson 
exchange potentials. 
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Figure 5: Real and imaginary part of the NN-scattering amplitude J{p, q) generated by iterating 
the (static) irreducible one-loop potential with the one-boson exchange. For further notations, 
see Fig. 2. 




Figure 6: Crossed three-boson exchange diagrams. The dashed lines represent the isovector- 
scalar boson of mass m. 



that modulo their isospin factors 7ti • T2 ± 6 the irreducible part of the right (3-rung ladder) 
diagram ^ in Fig. 1 agrees with that of the left diagram in Fig. 6 [8] . Discarding again the isospin 
factors 7ti ■ 72 ± 6 and — Ti ■ 72 ± 6 we deduce furthermore that the total sum of the three crossed 
three-boson exchange diagrams in Fig. 6 is equal to the difference between the right (middle) 
crossed diagram in Fig. 6 and the irreducible part of the right (middle) ladder diagram in Fig. 4. 
Putting back the isospin factors and using the four relations between the irreducible parts of the 
six two-loop diagrams, we find for the isoscalar component of the (static) three-boson exchange 
potential in momentum space: 

Vaiqr" = -g" / , 2 4 2 = arctan — , 13 

where we have made use of the results derived in the appendix of ref.[8]. The quantities loj — 



-I- Ij denote the on-shell energies of the three exchanged scalar-isovector bosons, where 

^3 = h + I2 + q- The corresponding potential in coordinate space is repulsive and can be 
expressed by a simple exponential function: Vs{rY'^^^ — 3g'^e~^"*^/(327r^mr^). The isovector part 
of the (irreducible) three-boson exchange potential proportional to ti • T2 reads on the other 
hand: 



(14) 
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At zero momentum transfer q = the six-dimensional integral in Eq.(14) can be evalu- 
ated analytically and the result (67rm)"~^ is found. For a numerical evaluation of the isovector 
potential V3(g)*^*''^ the representation in Eq.(14) is not best suited. A better representation is 
obtained by introducing Feynman parameters for each four- dimensional loop integral and per- 
forming the elementary integrations over some of the Feynman parameters. After several skilful 
transformations we arrive at the following very handy double-integral representation: 
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In this form the dx-integral could even be solved in terms of square-root and arctangent or 
logarithmic functions (depending on the sign of a radicand). The associated coordinate space 

potential Vy,{rY'^'"'> can now also be computed easily. One just has to replace the polynomial inside 
the square brackets in Eq.(15) by the expression —4:'Krx{l — x)ycx^{mr[{x'^ + xy + y + y'^)/{x{l — 
a;)|/)]^/^}. In the first and second row of Table 1 we have collected numerical values which display 
the dropping of the isoscalar and isovector three-boson exchange potentials 1^3(5)^**''^^ with the 
momentum transfer q. In each case we have divided 1^3(0')^**''"^ by the corresponding value at 
g = 0. For comparison, we give in the third row the analogous ratios for the one-loop two- 
boson exchange potential V2{q) written down in Eq.(9). One observes that the isovector three- 
boson exchange potential drops somewhat faster than its isoscalar counterpart. This decrease 
is however weak in comparison to that of the two-boson exchange potential, mainly due to the 
different mass scales (3m versus 2m) involved. As an aside, we note that for a scalar-isoscalar 
boson the (static) irreducible two-boson and three-boson exchange potentials would both vanish 
identically: V2{q) = V^lq) = 0. The latter zero-result follows from the relations between the 
irreducible parts of the six three-boson exchange two-loop diagrams, mentioned above. 
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3-isoscalar 


0.965 


0.882 


0.785 


0.695 


0.618 


0.554 
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0.416 


0.384 


3-isovector 


0.957 


0.856 


0.739 


0.633 


0.544 


0.471 


0.412 


0.363 


0.322 


0.289 


one-loop 


0.861 


0.623 


0.442 


0.323 


0.245 


0.192 


0.154 


0.127 


0.107 


0.091 


2-isoscalar 


0.927 


0.785 


0.655 


0.554 


0.476 


0.416 


0.369 


0.331 


0.300 


0.275 


2-isovector 


0.960 


0.874 


0.783 


0.704 


0.637 


0.581 


0.534 


0.494 


0.460 


0.430 



Table I: First and second row: The isoscalar and isovector three-boson exchange potentials 
Eqs.(13,15) versus the momentum transfer q. The respective values at q — have been di- 
vided out. The third row gives the analogous ratios for the one-loop two-boson exchange potential 
Eq.(9). The fourth and fifth row correspond to the two-boson-exchange potentials with one-loop 
vertex corrections written in Eqs.(18,19). 




Figure 7: Two-boson exchange diagrams with one-loop vertex corrections. The diagrams with 
uncrossed boson lines and the role of both nucleons interchanged are not shown. 



For the sake of completeness (at two-loop order in the nonrelativistic approximation) we 
consider also the two-boson exchange diagrams with vertex corrections shown in Fig. 7. In order 
to evaluate them we make use of the methods outlined in ref.[9]. From the one- loop corrections 
(beyond mass and coupling constant renormalization) to the isospin-even and isospin-odd boson- 
nucleon scattering amplitudes:^ 
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one can calculate (via unitarity) the spectral functions Im V2(i/i). These mass spectra determine 
then the two-boson exchange potentials in momentum space through an unsubtracted dispersion 
relation: 
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We find for the isoscalar part of the two-boson exchange potential with one-loop vertex correc- 
tions: 
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whose associated coordinate space representation is a simple exponential function: V2(r)['c^ = 

^The systematic expansion in powers of 1/M gives rise also to a (complex- valued) piece of order which is 

to be interpreted a relativistic l/M^-correction to the twicc-itcratcd boson-exchange. 
^For a scalar-isoscalar boson these one-loop corrections would vanish identically. 



o 



_3g,6 g 2mr j (iSTT^mr^) . The isovector component proportional to Ti • reads on the other hand: 



with the value —g^/ (27r^m^) at zero momentum transfer g = 0. In order to obtain the associated 
coordinate space potential one just has to replace arctan(g/2mx) in Eq.(19) by the expres- 
sion — e~2"^''^/ (47rr^). The numbers in the fourth and fifth row of Table I display the dependence 
of these two-boson exchange potentials with one-loop vertex corrections V2{q)^c'^^^ 
mentum transfer q. One observes a slow decrease similar to that of the irreducible three-boson 
exchange potentials. 

In summary we have calculated in this work all nonrelativistic contributions from scalar- 
isovector boson-exchange between nucleons at one- and two-loop order (including vertex cor- 
rections). At one loop order these contributions consist of the once-iterated boson-exchange 
amplitude and the irreducible two-boson exchange potential. At two-loop order one has the 
complex-valued scattering amplitudes related to the twice-iterated boson-exchange and the iter- 
ation of the one-loop potential with the one-boson exchange as well as the irreducible three-boson 
exchange potentials and the two-boson exchange potentials with one-loop vertex corrections. In 
each case we have given analytical expressions involving only the minimum number of integra- 
tions (to be done numerically). The present results could be useful for phase-shift analyses, few- 
and many-body calculations, etc. The applied methods can be straightforwardly generalized to 
the pseudoscalar-isovcctor pion with its spin- and momentum-dependent couplings to the nu- 
cleon. The number of amplitudes and diagrammatic contributions to each amplitude increases 
however drastically by several orders of magnitude. Work along these lines is in progress [7]. In 
passing we have to point out that there exist still other two-loop contributions which also scale 
as g^M^. These complex-valued amplitudes correspond either to the iteration of the relativistic 
1/M-correction to the one-loop potential with the one-boson exchange, or to relativistic 
correction to the twice-iterated boson-exchange. More work is necessary in order to isolate, 
classify, and analytically compute these relativistic correction terms. 
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